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Abstract. We develop a detailed model of the Milky Way (a "prototypical' disk 
galaxy) and extend it to other disks with the help of some simple scaling relations, 
obtained in the framework of Cold Dark Matter models. This phenomenological 
("hybrid") approach to the study of disk galaxy evolution allows us to reproduce 
successfully a large number of observed properties of disk galaxies in the local 
Universe and up to redshift z ~1. The important conclusion is that, on average, 
massive disks have formed the bulk of their stars earlier than their lower mass 
counterparts: the "star formation hierarchy" has been apparently opposite to the 
"dark matter assembly" hierarchy. It is not yet clear whether "feedback" (as used 
in semi-analytical models of galaxy evolution) can explain that discrepancy. 
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1. Introduction 

In the currently popular scenario for galaxy formation, pioneered by 
White and Rees (1978), the dark haloes of galaxies form hierarchically 
by the gravitational clustering of non-dissipative dark matter, while the 
luminous parts form through a combination of gravitational clustering 
and dissipative collapse (which may be affected by feedback). The ma- 
jor uncertainty in this scenario concerns the baryonic component, since 
the physics of star formation and feedback are very poorly understood 
at present, thus requiring a parametric approach to the problem. Semi- 
analytic models of galaxy formation help to explore large regions of 
the relevant parameter space and have produced quite encouraging 
results (e.g. Cole et al. 2000). However, they have not yet managed 
to reproduce succesfuUy some key observed properties of disks (van 
den Bosch 2002, see Sec. 4) or ellipticals (Peebles 2002, Thomas et 
al. 2002) [Note: see also the discussion report by Matteucci in these 
Proceedings]. 

Simple (i.e. not dynamical) models of (chemical and/or photomet- 
ric) galaxy evolution use, in general, fewer free parameters than semi- 
analytical models. If the number of observables that are successfully 
reproduced is (much) larger than the number of free parameters, it 
is reasonable to assume that the galaxian histories produced by such 
models may indeed match the real ones. In that spirit, we have de- 
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Figure 1. Chemical {left) and photometric {right) evolution of the Milky Way disk, 
according to the model of Boissier and Prantzos (1999). In all panels the solid curves 
correspond to model results at galactic ages of 1, 4 and 13 Gyr, respectively; the 
latter {heavy curves) are compared to observations of the present day disk (in the 
left panels: shaded regions for the gaseous and stellar profiles and data points for the 
Star Formation Rate). The model leads naturally to different scalelengths for the 
B-band (4 kpc) and the K-band (2.6 kpc), in agreement with observations. 

veloped a simple approach to disk galaxy evolution, based on i) a 
detailed model of the Milky Way (used as a prototype, Sec. 2) and ii) 
an extension to other disks trough some simple scaling relations (Sec. 
3). The main conclusion (Sec. 4) is that massive disks have formed 
the bulk of their stars earlier than low mass ones. Whether this "star 
formation hierarchy" is compatible with the "dark matter assembly" 
hierarchy remains to be demonstrated. 



2. A simple model for the Milky Way disk 

Several simple models have been developed in the past few years on the 
chemical evolution of the Milky Way disk (e.g. Prantzos and Aubert 
1995, Chiappini et al. 1997, Boissier and Prantzos 1999, Chang et al. 
1999, Portinari and Chiosi 2000). Although they may differ consid- 
erably in their ingredients (radial dependence of Star formation rate 
and/or infall rate, inclusion or not of radial inflows, etc.) these models 
converge in the following points (see e.g. discussion in Tosi 2000): i) 
infall (of primordial or low metallicity gas) on a long timescale (~7 
Gyr) is required locally, in order to reproduce the observed G-dwarf 
metallicity distribution in the solar neighborhood; ii) strong radial de- 
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Figure 2. Results of the Milky Way disk model of Boissier and Prantzos (1999) 
plotted as a function of the redshift (for a Universe with JIa/^O.SO, nA=0.70 and 
Ho=65 km/s/Mpc). The model suggests that such disks are fainter in the K-band 
in the past, but slightly brighter in the B-band at intermediate redshifts; they are 
also more compact in the past, especially in the B-band. 



pendence of the Star formation rate and/or the infall rate is needed in 
order to reproduce the observed radial profiles of gas, SFR and metal 
abundances [Note: in recent cosmological hydro-simulations of galaxy 
formation Sommer-Larsen et al. (2002) find that gas accretion in regions 
8 kpc distant from the centers of Milky Way type spirals takes place on 
timescales of ~6 Gyr, in reasonable agreement with the requirements 
of simple models] . 

Fig. 1 displays the results of a simple model, (from Boissier and 
Prantzos 1999) which also reproduces all the main observables in the 
solar neighborhood. The adopted SFR as a function of radius R is of 
the form 

^{R) = A V{R) R-^ ^gHR) (1) 

where V{R) is the rotational velocity at R and the gas surface 
density; this SFR is based on the assumption that large scale star 
formation in spirals is induced by the passage of the spiral waves (Wyse 
and Silk 1989). 

The model reproduces satisfactorily the observed chemical and pho- 
tometric radial profiles of the MW disk (Fig. 1) as well as those of 
oxygen and other metal abundances (Hou et al. 2000); in particular, the 
model predicts a flatenning of the abundance gradients with time (due 
to the assumed inside-out formation of the disk) , in agreement with the 
recent observational results of Maciel and da Costa(2002), based on O 
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abundances of planetary nebulae of various age classes [Note: taking 
into account the uncertainties in evaluating ages of planetary nebulae, 
the importance of that agreement should not be overestimated; it is, 
however, encouraging]. 

Fig. 2 displays the evolution of several observables of a Milky Way 
disk as a function of redshift in a "standard" cosmology, according 
to the model. It appears that such disks evolved with quasi-constant 
amount of gas and were in the past substantially more compact in 
the B-band; they also were slightly brighter in the B-band (at least at 
redshifts z ~ 1-2) and always fainter in the K-band. 

Although the model reproduces successfully the present day profiles 
of the Milky Way disk, its "predictions" for the past history of such 
disks cannnot be considered very robust, since a degenerate solution 
(several different histories leading to the same final state) cannot be 
excluded. Still, it is tempting to explore the consequences of that model, 
after extending it to other disk galaxies. 



3. Scaling relations for galactic disks 

Assuming that the MW is a typical spiral, one may try to extend the 
model to the case of other disk galaxies by means of simple scaling 
relations (in an analoguous way, the Sun is a typical star used to 
"calibrate" models of stellar evolution). The scaling relations derived 
by Mo et al. (1998) in the framework of Cold Dark Matter models of 
galaxy formation can be used (Boissier and Prantzos 2000, BP2000) 
to describe disks as two-parameter family (scale length Rd and central 
surface brightness Sq), in terms of the corresponding MW parameters 
(represented by subscript G below): 

Rd A Vc 



and 



RdG Ag Vcg 

md / A \ Vc 



(2) 
(3) 

SoG rndG \>'G/ Vcg 

The two fundamental parameters entering these relations are the circu- 
lar velocity Vc (a measure of the disk's mass) and the spin parameter 
A (a measure of the disk's angular momentum distribution); large Vc 
values correspond to massive disks and large A values to extended 
ones. The baryonic fraction rud of the disks is taken to be 0.05 in our 
calculations. Finally, for the MW parameters we have: i?dG=2.6 kpc, 
EoG=1150 Mq/pc^, VbG=220 km/s. 
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Note that these scahng relations apply to gaseous disks, assumed to 
be formed at the redshift of the corresponding dark halo formation; this 
formalism is adopted in semi-analytic models of galaxy evolution, where 
the baryon accretion history (the "infall" history, in our terminology) is 
determined essentially by the dark matter accretion history. A Schmidt 
type law for star formation (with an efficiency proportional to local 
dynamical timescale) is usually assumed in such models, along with 
some prescription for "feedback" (a term describing the way gas is 
affected by the energy input from stars) ; the interplay between the three 
effects (accretion, star formation, feedback) determines the effective 
star formation efficiency in semi-analytic models. 

A different approach is adopted in BP2000: the star formation law 
and efficiency A (Eq. 1) are kept the same as in the MW model (with 
rotational velocity curves properly calculated), while the infall rate is 
adjusted in order to reproduce main properties of present-day spirals. 
This approach (i.e. using present-day properties of galaxies to infer 
their histories) has been characterised as "backwards" approach to 
galaxy evolution (e.g. Fereras and Silk 2001); our use of scaling relations 
borrowed from Cold Dark Matter ("forwards") models qualifies our 
approach rather as a "hybrid" one. In fact, the scaling relations allow 
to tie the geometric properties of disks to those of the Milky Way and 
nothing more; the overall disk evolution is really determined by the 
assumed infall (or mass accretion) history. 



4. Results 

Our "hybrid" approach allows us to reproduce successfully a large num- 
ber of observables of galaxies in the local Universe (BP2000), including: 
disk sizes and central surface brightness, Tully-Fisher relations in vari- 
ous wavelength bands, colour-colour and colour-magnitude relations, 
gas fractions vs magnitude, abundances as a function of local and 
integrated properties etc; it also reproduces naturally detailed spectra 
from the Kennicut atlas of galaxies for various galaxy types, as can be 
seen in Fig. 3. 

We find that the crucial ingredient of that success is the assumption 
that gas is infalling more slowly in galaxies with smaller masses and/or 
surface densities. These two parameters are the main drivers of disk 
evolution. We note that Bell and de Jong (2000) and Bell and Bower 
(2000) reach similar conclusions, but they attribute a more important 
role to local surface density while we find that galaxy mass plays an 
even more important role: on average, massive disks are older than 
lower mass ones. Disks with Vc ~100 km/s have formed the first half 
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Figure 3. Synthetic spectra {solid curves, absorption features only) of model disks 
at an age of 13 Gyr compared to observed spectra of various Hubble types {dotted 
curves, with emission features) from the Kennicutt atlas. Model spectra are given for 
low (upper panels), intermediate (middle panels) and high (bottom panels) circular 
velocities, and for low (left panels), intermediate (middle panels) and high (right 
panels) values of the spin parameter A. Each model galaxy is compared to an 
observed one of approximately equal circular velocity Vc- It can be seen that more 
massive disks are, in general, redder, but also that for similar values of Vc the 
value of A (i.e. surface density) is important for the age (and colour) of the stellar 
population (from Boissier and Prantzos 2000). 

of their stars within 4-9 Gyr (Fig. 4), while disks with Vc ~300 km/s 
have done so within 2-6 Gyr (depending on spin parameter, that is 
surface density, see Boissier et al. 2001). 

The main observational arguments on which our conclusion is based 
are presented in Fig. 5. One sees that, on average, massive disks are 
redder, poorer in gas and richer in metals than less massive ones. Each 
one of these observables, taken separately, could be easily explained 
(e.g. redder colours could be explained by larger amounts of dust, 
as invoked e.g. in Somerville and Primack 1999 or Avila-Reese and 
Firmani 2000). When all three observables are considered together, the 
idea of a star formation efficiency increasing with galaxy mass apears as 
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Figure 4- Time required for formation of the first 50% of the stars of a disk, according 
to models of Boissier and Prantzos (2000), as a function of circular velocity. The 
three curves correspond to spin parameters A equal to 1/3, 1 and 3 times the one of 
the Milky Way disk. Small mass and/or extended disks form their stars on longer 
timescales, i.e. both mass and surface density affect the formation timescale. 

a viable explanation. Such a possibility has been invoked in Ferreras and 
Silk (2001), but it requires huge amounts of (undetectable) gas to be 
present in small galaxies. However, when the fourth observable is taken 
into account, namely the absence of any dependence of the present-day 
star formation efficiency on rotational velocity (or B-magnitude, see 
Boissier et al. 2001), the only viable explanation is the one invoked 
above: massive disks are older, on average. 

We note that in recent cosmological hydro-simulations Nagamine et 
al. (2001) find that star formation in small galaxies has stopped many 
Gyr ago, in clear contradiction with observations of local galaxies. On 
the other hand, van den Bosch (2002) finds that semi-analytical mod- 
els, even with feedback, produce massive disks that are systematically 
bluer than their lower mass counterparts, again in contradiction with 
observations (Fig. 6); the reason of the failure is obviously related to the 
fact that the mass accretion histories of baryons are largely dictated by 
the hierarchical clustering of dark matter (e.g. Avila-Reese and Firmani 
2000). Once gas becomes available it forms rapidly stars; feedback can 
only delay star formation for a short time (shorter than the several 
Gyr that are observationally required to obtain small disks bluer than 
massive ones). 

Our simple, "hybrid" model for disk evolution, calibrated on the 
MW, suggests that on average massive disks have formed the bulk of 
their stars several Gyr earlier than low mass ones. Their predictions 
match successfully most currently available observables, including data 
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Figure 5. From top to bottom: Colour, metallicity at 3 kpc from the center, gas 
fraction and star formation efficiency (defined as star formation rate divided by the 
total gas amount) as a function of rotational velocity for galactic disks. In general, 
more massive (rapidly rotating) disks are redder, more metal-rich and gas poor than 
less massive ones; however, the star formation efflciency does not seem to depend 
on the disk's mass. Put together, these observational data suggest that the more 
massive disks are chemically and photometrically older than less massive ones. The 
curves correspond to model results of Boissier et al. (2001) for three values of the spin 
parameter A: 1/3, 1 and 3 times the one of the Milky Way (the latter corresponds 
to the solid curves in all panels). 



from surveys at intermediate redshifts: as shown in Boissier and Prant- 
zos (2001) this simple model may account for the lack of evolution in 
the properties of large disks observed up to z ~1 by the Canada- France- 
Hawaii survey (Lilly et al. 1998): the evolution of large (and massive) 
disks has been achieved mostly before z ~1 in our senario. 

A prediction of our model appears in Fig. 8. The anticorelation 
between the metal abundance gradient and the inverse B-scalelength, 
found to be valid locally (Prantzos and Boissier 2000)is shown to be 
valid also at high redshifts: smaller disks are always caracterised by 
larger (negative) abundance gradients. 
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Figure 6. Colour vs. Magnitude relation of galaxies obtained by semi-analytic mod- 
els without feedback (left) and with feedback (right) by van den Bosch (2002). The 
results (points) are compared to the observational trend (straight lines) ; there is an 
obvious disagreement of models with observations. 



5. Conclusion 



In the currently popular paradigm of hierarchical galaxy formation, 
low mass dark matter haloes form first, while more massive ones are 
formed later through accretion and merging; in principle, baryons are 
supposed to follow the dark matter, but their fate is largely unknown 
at present, due to a lack of a reliable theory of star formation (and 
feedback) . 
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Figure 7. Properties of large (B-scalelength >4 kpc) disks as a function of redshift. 
Data points are from the CFH survey (Lilly et al. 1998). The shaded aereas corre- 
spond to our models (1 ct around the mean value for the dark shaded aerea and 3 
a for the light shaded one), with the solid curves representing the mean value and 
the dotted one the Milky Way evolution. 
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Figure 8. Oxygen abundance gradient vs the inverse of the B-band scalength in disks 
for 3 different values of redshift. The anticorrelation found for local disks (shaded 
aerea on the left panel) by Prantzos and Boissier (2000) and which is supported by 
observations (data points on the left panel), is also found at higher redshifts. 



At present, and despite claims to the contrary, there is no satis- 
factory explanation (at least, not a published one) for the observables 
presented in Fig. 5 in the framework of hierarchical galaxy formation. It 
remains to be shown why star formation in galaxies apparently followed 
an "inverted hierarchy" w.r.t the dark matter asembly. Feedback offers 
an obvious solution to that problem, but the required delay timescales 
appear unphysically large. On the other hand, we have shown that 
simple models based on our (fairly detailed) knowledge of the Milky 
Way, reproduce most properties of local disks with few parameters 
(essentially the infall timescale) and make interesting predictions for 
the disk properties at higher redshifts. 
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